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We discuss spin-^ two-dimensional (2D) Heisenberg antiferromagnet (AF) on a square lattice 
at T = in strong magnetic field H near its saturation value He. A perturbation approach is 
proposed to obtain spectrum of magnons with momenta not very close to AF vector in the leading 
order in small parameter (He — H) / He. We find that magnons are well-defined quasi-particles at 
H > 0.9He although the damping is quite large near the zone boundary. A characteristic rotonlike 
local minimum in the spectrum is observed at k = (vr, 0) accompanied by decrease of the damping 
near (7r,0). The suggested approach can be used in discussion of short- wavelength excitations in 
other 2D Bose gases of particles or quasi-particles. 

PACS numbers: 75.10.Jm, 75.50.Ee, 75.40.Gb 

I. INTRODUCTION 

Spin-^ two-dimensional (2D) Heisenberg antiferromagnet (AF) on a square lattice has been one of the most attrac- 
tive theoretical objects in the last two decades because this model describes parent compounds of high- Tc supercon- 
ducting cuprates^. Whereas the theory of long- wavelength magnons in quantum square 2D AF is well developed and 
describes well existing experimental data^^^^^, there are some surprising recent experimental findings indicating that 
the standard theoretical approaches do not work for short- wavelength magnons. Thus, a rotonlike local minimum was 
observed at small T in the spin- wave dispersion at k = (tt, 0) in a number of recent experiments on square spin-^ 2D 
^Pg3j4j5^,6 ^]^ig minimum is not reproduced quantitatively within second order of 1/S expansion and phase flux RVB 
techniques, while numerical computations using Quantum Monte-Carlo and series expansion describe it satisfactorily 
(see Ref.^ and references therein). The origin of the rotonlike minimum has not been clarified yet. It is attributed to 
the entanglement of spins on neighboring sites^. Rotonlike minimum in the spectrum of short- wavelength magnons 
(at k = (tt, 7r/>/3)) has been reported recently also in spin-^ triangular AF at i7 = obtained using 1/S expansion^i^ 
and series expansion^' 

Even more peculiar finding concerning short- wavelength magnons in square spin-^ 2D AF was obtained in Ref.— . 
The authors investigated renormalization of the spin- wave spectrum within first order in 1/S in strong magnetic 
field H smaller than the saturation field He. Due to noncollinearity of sublattices and the field, three-magnon terms 
appear in the Hamiltonian after transformation of spin operators to bosons. Such terms lead to the loop diagram 
for the self-energy part S(/c) in the first order in 1/S corresponding to processes of magnon decay into two magnons. 
The spin- wave damping obtained within the first order in 1/S was not very large. To discuss decaying processes 
self-consistently authors went beyond the first order in 1/S carrying out a self-consistent calculation of the Green's 
function (GF) G(a;,k) using the dressed GF considering the loop diagram. As a result they obtained that ImG(a;,k) 
as a function of oo at fixed k does not resemble a peak at all in most of the Brillouin zone (BZ) when 0.76Hc < H < He 
(the worse situation was dit H ^ 0.9 He). Thus the authors concluded that magnons are unstable according to decay 
into two magnons at H ^ He. Magnons reappeared only at > 0.99 He due to decreasing of the three-magnon 
vertex. Meantime within the approach suggested in Ref.^^ only some higher-order 1/S corrections are taken into 
account and it is not clear what is the role of the ignored 1/S terms. This question looks very serious because there is 
no small parameter in 1/S expansion for S = 1/2. On the other hand the unsuccessful attempts to describe rotonlike 
minimum in the spectrum within second order in 1/5 at = and results of Ref.^^ themselves signify that higher 
order 1/S terms can renormalize the spectrum of short- wavelength magnons noticeably if 5 ~ 1. 

In view of the above mentioned results it is highly desirable to develop an analytical perturbation theory allowing 
quantitative discussion of the short- wavelength magnons in quantum 2D AF. The aim of the present paper is to 
develop such theory in the parameter {He — H)/He <C 1 for magnons with momenta not very close to AF vector 
ko = (7r,7r). The magnon spectrum is calculated below in the leading order in {He — H)/He. Our results can be 
tested experimentally on the growing family of AF compounds with small or moderate exchange coupling constant 
J (i.e., with experimentally achievable He cx J) which can be modeled by square spin--^ Heisenberg AF^l^^. The 
approach suggested below can be used in discussion of short- wavelength excitations in other 2D Bose gases of particles 
or quasi-particles. 

The rest of the paper is organized as follows. We discuss the approach in Sec. [Ill Results of the spectrum calculation 
are presented in Sec. Illli Sec. IIVI contains our conclusion. One appendix is added devoted to diagrams estimations. 
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FIG. 1: (Color online.) (a) Equation for the vertex r(p, /c), where k = (ct^,k). (b) Expression for the self-energy part in 
the leading order in 11/ He for k not very close to ko and (7r,0). (c) and (d) Some higher order diagrams for which are 

much smaller than those presented on slide (a) ii /i <^ He and k is not very close to ko and (7r,0). Lines with two arrows 
and wavy lines correspond to anomalous Green's functions and condensed particles, respectively. Diagrams the more simple of 
which are shown on slide (d) contribute to the leading correction to ReS(/c) at k (tt, 0). 



II. MODEL AND TECHNIQUE 



We discuss spin-^ Heisenberg AF on the square lattice in magnetic field which Hamiltonian has the form 

W = ^E'^^iSiS, + ff5]5f. (1) 

We express spin projections via Pauli operators a\ and ai as Sf = — ^ +ajai, S\ = aj, 5*^" = a^. To treat them as 
Bose operators one should introduce to Hamiltonian ([T]) the constraint term U a\a\aiai^ where ^ 00, describing 
infinite repulsion of particles on the same site^^. As a result Hamiltonian ^ reads as 

H ^ (ek - /i) aj^ttk 

k 

+ ^ (^ki-k3 +^)4i42^k3«k4, (2) 

ki+k2=k3+k4 

where ek = ^(^0 + ^k), ^k = Xli Jij cos(kR^j), R^j connects sites i and j, 11 = He — H ^ He = Jo is the saturation 
field, and N is the number of spins in the lattice. The bare spectrum ek is quadratic near ko. 

Formally Eq. ([2]) describes a dilute Bose gas with pare interaction between particles and with the chemical potential 
ji. In 3D AF this equivalence leads to the solution at /i <C Hc^^ by borrowing the results from the theory of the dilute 
Bose gas^^'^^. Within this approach the long-range AF ordering in the plane perpendicular to the field appearing at 
H < He corresponds to condensation of magnons with momentum equal to ko. The density of particles is proportional 
to jii = He — H . Therefore one can use results for dilute Bose gas when /i <C i^c- The ladder approximation for the 
vertex is valid and one has for the vertex the equation shown in Fig. [TJa). Normal and anomalous TL{k) self- 
energy parts are expressed via the vertex, where k = (c^,k) and the same notation is used below. In particular one 
has dit H < He for the expression shown in Fig.llfb) in the leading order in fi/He and corrections to it are small 
from other diagrams some of which are shown in Figs.[l](c) and[Hd). 

It is tempting to apply the same approach to 2D AF at T = but in contrast to 3D dilute Bose gas the theory of 
2D one at T = can be developed only in the limit of exponentially small density of particles^^. It means in the case 
of 2D AF that there is a good perturbation theory only for fields exponentially close to i^c, ln{He//J^) ^ 1, while it 
is desirable to extend the theory to a more acceptable smallness of the parameter 11/ He. We show now that one can 
overcome this obstacle and construct a perturbation theory on fi/He to obtain the magnon spectrum for momenta 
|k — ko| ~ 1 in the leading order. 
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FIG. 2: (Color online.) A quarter of Brillouin zone is shown and an area around ko is drawn in red in which the suggested 
perturbation approach does not work when hi{Hc/ fj) ^ 1 (i.e., in which the condition (|3]) does not hold). The red area grows 
as the field reduces. 



We remind first that self-energy parts are small at /i <C i^c and they vanish at /i = Oi2.dI. Then, one notices that 
for such k that the inequality 

ek-/i> |S(ek,k)| (3) 

satisfies, one can neglect |n(/c)p in the GF denominator T){k) = (cj — ek + /i — Ti{k)){uo + Ck — /i + S(— /c)) + |n(/c)p 
and the normal GF G{k) = (cj + ek - /i + T.{-k))/V{k) (see RefsJ^^ii) has the form 

It is implied in Eq. (|3]) that |i;(ek,k)| ~ |n(ek,k)| that is really the case in the entire BZ as particular estimations 
show. Inequality ([3]) is not satisfied only in an area around ko that grows as the field reduces (see Fig. [2]). Thus one 
has to calculate only normal self-energy part to obtain the magnon spectrum at |k — ko| ~ 1 (to be precise, when the 
condition (j3|) holds). 

It is the matter of direct estimations, similar to those carried out in RefsJ^'— for dilute 3D Bose gas, to show using 
the smallness of S and H that we have expression for S(/c) presented in Fig. [IJb) in the leading order in ji/ He if k 
is not very close to ko and (tt, 0) (see Appendix for some details). Integrals in other diagrams, some of which are 
presented in Figs.fljc) andfT^d), become less singular at small ji at such external momenta as compared with cases of 
k ~ ko and (7r,0). As a result contribution from them estimated in the Appendix turns out to be of the next order 
in ji/ He. Diagrams, the more simple of which are shown in Fig.[T](d), should be also taken into account at k ~ (tt, 0) 
as we discuss below. 

To calculate directly two diagrams drawn in Fig. [T](b) one needs to know the condensed particle density po and 
the spectrum in the entire BZ which can be found only using self-energy parts near ko^^^^^. The trick here bases on 
the fact that momenta of integration near ko are essential in the second diagram in Fig. [Tfb). Thus, one can put the 
momentum q = {uq, q) in the vertex equal to ko = (0, ko) and factor out the vertex. As a result we have 

S(/c) = 4nr(0,/c), (5) 

where t ^ -\-0 and n is the density of particles which is expressed in spin-^ AF via the magnetization M = — "^^{S^) 
as n = ^ — The problem remains to calculate M but we can make use of results of previous numerical calculations. 

Magnetization in spin-| square 2D AF in magnetic field was investigated before in the second order in and 
by numerical diagonalization of finite clusters^. It was obtained in Ref.^^ that expression for the magnetization 
derived within the first order in 1/5 fits the numerical data very good except for close vicinity of He. Thus we can 
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approximate n to high accuracy by its value obtained in the first order in 1/S for which one has after calculations 
following those of Ref 



/i ^ 1 - /i Jq / Jq + Jq 



2 TV ^ 2 Jo y Jo + Jq - 2 Jq/i(2 - fl) ' 



(7) 



q 



where 



where ji = ji/ He. In particular, one obtains from Eq. ([7j) n ~ 0.083 and 0.048 for H = 0.9 He and 0.95i7c, respectively. 
We conclude from Eq. ([7j) that n ^ fi/ He if H is not exponentially close to He. 

It remains only to obtain the vertex r(0, /c). We imply for simplicity that there is exchange coupling between 
neighboring spins only with the exchange coupling constant J. It is convenient to try solution of the equation shown 
in Fig. [TJa) in the form 

r(p, k) = a{k) + P{k) sin (^p^ + ^k^^ + -f{k) sin (^p, + ^k,^ . (8) 

Putting this expression in the equation one leads to a set of three linear algebraic equations on /3{k) and ^{k) 

that can be readily solved. The result is quit cumbersome and we do not present it here. One leads to very compact 
expressions for the vertex at p = and k = {u = -\- iS, k) in the particular cases of \kx\ = \kz\ and kz = having 
the form 

T{0,k) = for \k,\ = \U (9) 

'^(k+ko)/2 "^0 

T=-T ^ . (11) 

TV ^ 2(6q + eq+k+ko - ek - 

We use these results below to calculate for k directed along vectors kg and (7r,0). As we discuss corrections to 
the spectrum in the leading order in /i/iJc, we discard fi in the expressions for the vertex. 

It is difficult to calculate the vertex analytically in the entire BZ and we perform in the next section numerical 
integration to find it along some particular directions in BZ and for some particular H values. 

III. RESULTS OF THE SPECTRUM CALCULATION 

Results for the momentum scan along ko obtained using Eqs. ([9]) and (pT]) are plotted in Fig. [3] for iJ = 0.9iJc 
and 0.95iJc- Solid (dashed) line in Fig. [3] is for the area in which inequality ([3|) does (does not) hold. Momenta kc 
at which the solid line turns into the dashed one have been found from the assumption |S(ek^,kc)| = 0.2(ek^ — /i) 
that gives kc ~ 0.47ko and 0.6ko for H = 0.9He and 0.95iJc, respectively. Insets in Fig. [3] show that the ratio of 
the imaginary and the real part of the spectrum does not exceed 0.17 upon increasing |k| up to |kc|. The spin- 
wave spectrum is also drawn in Fig. [3] obtained within the linear spin- wave theory (LSWT) and having the form 
^y/{Jo + Jk)(Jo + (2(iJ/iJc)^ — l)Jk). Notice that at k = the spectrum obtained within our approach is real and 
coincides with that derived in LSWT. This finding is in agreement with quite general arguments presented in Refi^ 
according which quantum fluctuations do not change the classical spectrum at k = in quantum 2D AF at iJ 7^ 0. 
Quantum fluctuations shift down the spectrum at finite k as the first 1/S corrections dc^^. 

Let us discuss the spectral function defined as —^lmG{iJ + iS^k). As is explained above, the vertex r(0, /c), where 
k = (cj,k), can be found exactly using representation ([8|) with the following result for \kx\ = \kz\ 

^ -{a; - Jo)J§{2Tu; + (1 - 2T) Jq) + M-4Tu; + (-3 + 4T) Jq) Jf^^^^^^, - 2TJl^^^^y, 

^ ' ' 2Jo(2Tu;Jo + (1 - 2T) J2 - 2Tj2^^^^)/2) ' 

where T = 1/N Jo/2(eq + eq+k+ko — u; — iS). These expressions can be simplified further at a; = ek and one leads 
to Eqs. ([9]) and (pT|) . We plot in Fig. [H the spectral function for some momenta k||ko obtained using Eq. (p!2|) . One 
can see that curves look like peaks while peaks widths are quite large near the zone boundary. It should be noted 
that Fig. m is in contrast to the corresponding figure of Ref.^^ for H = 0.9He obtained using 1/S expansion. 
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FIG. 3: (Color online.) Real and imaginary parts of magnon spectrum (in units of the exchange coupling constant J between 
nearest spins) for k = <Jko calculated in the leading order in fi/Hc for (a) H = 0.9ifc and (b) H = 0.95i/c. Insets show ratios 
of imaginary and real parts of the spectrum. Dashed lines are for area in which the suggested perturbation technique does not 
work (inequality (|3]) does not hold). 

Results of the scan along another path in BZ are shown in Fig. [5l The most remarkable feature of this scan is the 
rotonlike minimum at (tt, 0) similar to that obtained recently at = Q^diL^ . It should be noted that imaginary part 
of T diverges at k = (tt, 0) and the vertex vanishes because one has in Eq. (pT]) Cq + eq+k+ko — = 2J(1 + cosqx). 
Notice also that for k||(7r, 0) the damping is at least of one order smaller than the real part of the spectrum even for 
H = 0.9Hc. As is noticed above, extra diagrams should be taken into account at k ^ (tt, 0) the more simple of which 
are presented in Fig. [Ifd) and which sum cannot be calculated analytically. It is shown in the Appendix that they 
give zero contribution to lmT^{k) at k = (tt, 0) in the leading order and one concludes that there is a decrease of the 
damping near (7r,0). At the same time their real part contributes to the leading correction to Rel](/c) at k ~ ('^r, 0) 
being of the order of /i. Then the real part of the spectrum can deviate near (tt, 0) from curves drawn in Fig.[5l These 
diagrams are of different signs at k = (tt, 0) (e.g., the first two diagrams in Fig. [Ifd) are negative and the last one is 
positive) and we cannot prove rigorously that there is the rotonlike minimum in the real part of the spectrum if H 
is not very close to He. Meantime we demonstrate now using the theory of 2D dilute Bose ga^ that the rotonlike 
minimum does exist if \n{Hc/ /i) ^ 1. One estimates in the leading order ReS(/c) ~ —ji/ hi{Hc/ jn) + Jk,^ \n{Hc/ fi) (the 
first and the second term here stem from the first diagrams shown in Fig. [Ifd) and[T](b), respectively) at k = (tt, k.) 
and K, <C jii/Hc while Rei;(A: ~ ■^) ~ /i In (i^c//i) and ek ^ 2 J — Jk^ /2. It seems to us likely that the minimum 
remains also at smaller H. 

It is difficult to go beyond the first order in ii/Hc because, in particular, one needs to know self-energy parts near 
ko to calculate the next order diagrams. We can estimate them to be 0{ii^/ijl/Hc) from the following consideration 
that works near (tt, 0) only for Iml](/c). Using relation S(/co) — n(/co) = /i one obtains for the spectrum near AF vector 
\/ek(ek + 2Il{kQ)Y^'^^ . As a result diagrams, some of which are shown in Figs. [TJc) andlljd), are of the order of 
pQ^J JIi{kQ) or smaller (see Appendix for details). On the other hand one concludes from Eqs. ([6]) and ([7j), bearing in 
mind that the second term in Eq. (|6|) is positive and of the order of n(/co), that po, n(/co) ~ /i if is not exponentially 
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FIG. 5: (Color online.) Same as in Fig. [3]but along the path in BZ shown in insets. Rotonlike local minima are seen at (tt, 0) 
accompanied by decrease of the damping near (7r,0). 



close to He. 

Thus, one can expect from estimations by the order of magnitude that the approach discussed in the present paper 
works for H > 0.9 He. It should be noted that numerical evidences have appeared recently that short- wavelength 
magnons are unstable in some regions of BZ at < 0.9Hc.^^ 
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IV. CONCLUSION 



In conclusion, we discuss spin-^ 2D Heisenberg AF on a square lattice at T = in strong magnetic field H near its 
saturation value He. A perturbation approach is proposed to obtain spectrum of magnons with momenta not very 
close to AF vector ko = (7r,7r) in the leading order in small parameter 11/ He = {He — H)/Hc. It is shown that only 
normal self-energy part contributes to the spectrum renormalization for which we have two diagrams shown 

in Fig. [Tfb). The sum of these diagrams is proportional to the number of particles which is expressed via uniform 
magnetization investigated before numerically (see Eqs. ([5j) and (|6])). These diagrams are of the order of /i. Higher- 
order diagrams some of which are shown in Fig.fljc) andfT^d) are of the order of 0{ii^ 1^1 He) at |k — ko| ~ 1 except 
for the vicinity of the point k ~ (tt, 0), where diagrams the more simple of which are shown in Fig.fljd) contribute to 
the leading part of Rel](/c). We find that magnons are well-defined quasi-particles at > 0.9i^c (see Figs. [3] and H]). 
A characteristic rotonlike local minimum in the spectrum is observed at k = (tt, 0) accompanied by decrease of the 
damping (see Fig. [5]). 

The approach suggested in the present paper can be used in discussion of short- wavelength excitations in 2D Bose 
gases of particles or quasi-particles. In the latter case one needs to calculate numerically the uniform magnetization 
to find the particle density n in Eq. ([5j). We note also that similar approach basing on the hard-core bosons formalism 
was proposed in Ref.^^ for AFs with singlet ground state. 
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APPENDIX A: DIAGRAMS ESTIMATIONS 



This Appendix is devoted to estimation of diagrams. We demonstrate first that when k is not very close to ko and 
(tt, 0) the leading order corrections to the normal self-energy part S(/c) come from two diagrams shown in Fig. [Tfb) 
and other diagrams are small some of which are shown in Figs. [Tfc) andlljd). As it is explained in the main text, 
diagrams presented in Fig. [ijb) are of the order of Jpo ^ 

Let us discuss the first diagram in Fig. [Hd). Representing the normal GF as G{k) = {E(k) + cj)/(cj^ — e^), where 
E{k) = ek — /i + S(— /c) and £k is the renormalized spectrum, one estimates for the first diagram in Fig. [TJd) after 
integration over the frequency and putting k = (^k^k) 



^(^) - (^(^ - A; + A;o) + Sq-k+ko)(^q + gq-k+ko) ^^^^ _ ^ (Al) 



^ q gq ^^q-k+ko[(gq + gq-k+ko)^ - 

One has in Eq. (|Aip when q ^ ko 

E{q)-e^^ |n(A:o)P ^ |n(A:o)| 



{E{q) + eq)eq (n(fco) + ek + £q)£q 



(A2) 



where we use the relation^iii ^(^0) — n(/co) = /i and Sq = \/eq(eq + 2n(/co)). Then, one has when q ^ ko and 
K > U{ko)/Hc: gq ~ Cq + n(/co) and 

gq + gq-k+ko - ^k ~ (/^^ - {f^l COS kx + 1^1 COS kz) " 2{i^x siu kx + i^z sin kz)) + n(/co), (A3) 

where k = q — ko- One concludes from the above consideration bearing in mind that E,n < /i that momenta of 
summation q near kg are essential in Eq. ()Aip . 

Let us discuss the point k = (|^, One has from Eq. (|A3p in this case ^q + Sq-k+ko — ^k ^ —J{nx + z^^) + Ii{ko). 
As a result one obtains from Eqs. (jAip . (|A2p . and (|A3p that the first diagram shown in Fig. [TJd) is proportional to 
iJ po \/ 11 1 He ~ iJ^iij Hcf*l^ ^ where we use also that n(/co) ^ as it is explained in the main text. 

One leads to the same estimations of Eq. ()Aip for other k except for k ~ (tt, 0) and k ~ ko- One has at k = (tt, 0) 
from Eq. ()A3p STq+Sq-k+ko^^k ^ JAi:^+n(/co). As a result the first diagram shown in Fig.[T][d) at k = (tt, 0) contributes 
to the main correction to Rel](/c) being of the order of J po ^ p. Its contribution to ImS(/c) is of the higher order in the 
parameter p/Hc. The same is true also for k ~ ko because we have at k = ko: s^q + s:q_k+ko ~^k = 2£q ^ J k!^ -\-Il{ko) . 
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We assume in the above estimations that r{k — ko^ko) ~ J in Eq. ()Aip . It is reahy the case as the results 
show of particular calculations of T{k — ko^ko) carried out as it is described in the main text. Meantime the vertex 
r(0,/c) <C J in the vicinity of the point k = (tt, 0) because T diverges (see the main text). Thus only diagrams 
containing T{k — ko^ko) contribute to the leading correction to ReS(/c) at k ~ (7r,0). The more simple of such 
diagrams are shown in Fig. [Tfd). 
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